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ABSTRACT 


In  this  note  an  exact  solution  is  given  to  the  problem  of  the  inter¬ 
section  of  an  arbitrarily  oriented  cone  and  a  sphere  for  the  application  of 
finding  the  Earth  coverage- -footprint- -of  a  narrow  beam  satellite  borne 
antenna.  Two  cases  are  considered:  when  the  pointing  angles  from  the 
satellite  are  known,  and  when  the  coordinates  of  the  point  on  the  Earth 
toward  which  the  boresight  is  aimed  are  known.  The  situation  of  partial 
coverage  is  also  treated  in  detail.  Examples  are  included  to  illustrate 
the  method. 
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List  of  Symbols 


a  =  vertex  angle  of  cone 
II  =  range  vector  from  satellite  to  Earth 
p  =  pitch  angle 
p  =  roll  angle 

_i,  j_,  k,  i_' ,  _I,  etc.  =  unit  vectors  in  respective  coordinate  systems 
6  =  instantaneous  declination  of  satellite 

r_  =  radius  vector  of  Earth 

C  =  curve  of  intersection  between  cone  and  sphere 
P  =  piercing  point  of  boresight  axis 
P  =  piercing  point  of  Mn  vector 
m  =  a  numeric 

=  longitude  separation  of  a  point  from  subsatellite  longitude 

u  =  an  auxiliary  vector  used  in  defining  surface  of  cone  and  vectors  Mn 

Mn  =  a  vector  on  surface  of  cone 

P  =  an  auxiliary  angular  parameter  used  in  locating  M 

—  n 

b,  c,  d  =  angles  of  triangle  formed  by  _i,  M_n  »  anc*  mr 

M  =  length  of  M  (M  ^  I M  I ) 
n  — n  n  1  — n 1 

T  =  tilt  angle  of  plane  of  triangle  relative  to  yz  -plane 
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EARTH  COVERAGE  (n FOOTPRINT n) 

OF  A  SATELLITE-BORNE  ANTENNA 

I.  INTRODUCTION 

The  purpose  of  this  Technical  Note  is  to  present  an  exact  solution  for 
the  problem  of  the  intersection  of  an  arbitrarily  oriented  cone  and  sphere. 

The  motivation  for  considering  this  problem  is  to  provide  an  aid  for  (a) 
evaluating  the  coverage- -footprint- -on  the  surface  of  the  Earth  of  a  satellite- 
borne  antenna,  and  (b)  assessing  the  threats  of  jamming  and  of  eavesdropping. 
Vectorial  and  trigonometric  methods  are  used  to  obtain  the  solution.  First, 
an  expression  is  derived  for  computing  the  latitude  and  longitude  coordinates 
of  the  point  on  the  Earth  where  the  boresight  axis  of  the  antenna  pierces  the 
surface  of  the  sphere.  Next,  a  chosen  number  of  vectors  are  generated  around 
this  axis  such  that  all  lie  on  the  surface  of  the  cone.  Another  expression  is 
derived  for  computing  the  coordinates  of  those  points  where  each  of  these 
vectors  pierces  the  surface  of  the  sphere.  A  computer  program  was  written 
for  solving  these  expressions.  The  outputs  of  the  program  are  a  listing  of 
the  coordinates  of  the  points  of  intersection  and  a  graphical  plot  of  these  points. 
Two  examples  are  included  to  illustrate  the  method. 

IIo  THE  INTERSECTION  OF  AN  ARBITRARILY  ORIENTED 

CONE  AND  A  SPHERE 

Consider  two  coordinate  systems:  OXYZ  and  Oxyz.  A  sphere  of 
radius  r  has  its  center  at  the  origin  of  the  first  one,  a  cone  of  vertex  angle 
a  has  its  vertex  at  the  origin  of  the  second  one.  The  orientation  of  one 
coordinate  system  relative  to  the  other  may  be  arbitrary;  however,  their 
instantaneous  relationship  is  known.  The  origins  are  separated  by  a  distance 
(m) (r )  measured  along  Oy,  where  m  is  a  numeric  and  it  is  greater  than 
unity  for  orbiting  satellites.  (The  simple  problem  of  intersection  of  a  cone 
and  a  sphere  when  the  vertex  of  the  cone  is  within  the  sphere  is  not  treated 
here.  )  The  range  vector  R  — which  is  also  the  axis  of  the  right  circular 
cone --from  the  origin  of  Oxyz  to  the  point  on  the  Earth  is  known.  Given  that 
a  cone  of  vertex  angle  a  is  described  around  axis  II  ,  one  would  like  to 
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know  the  curve  C  of  the  first  intersection  of  this  cone  with  the  sphere,  (The 
intersection  of  the  cone  emerging  from  the  sphere  is  of  no  interest.)  The 
problem  is  illustrated  in  Fig.  1. 

The  range  vector,  R_  ,  is  given  in  the  Oxyz  coordinate  system: 

R  =  A  R.  +  iRj  +  k  Rk  (1) 

where  _i»  J_»  k  are  unit  vectors.  Define  pitch  and  roll  angles,  p  and  p, 
as  shown  in  Fig.  1  : 

R. 

l 

sin  p  =  - 

1  J 

RK 

tan  p  = - - 

Vr.2  +  R.2' 

1  J 

The  R_  vector  can  be  written  in  terms  of  these  angles: 

II  =  _i  ( |  R.|  cos  p  sin  p)  +  ( |  R_|  cos  p  cos  p)  +  k  (  |  R_|  sin  p) 

(4) 

where  |  Rj  is  the  magnitude  of  R_  .  Let  Oxyz  be  rotated  around  Ox,  in 
the  direction  indicated  by  the  arrow  in  Fig.  1,  by  the  angle  6  ,  the  declina¬ 
tion  angle  of  the  satellite,  resulting  in  the  primed  coordinate  system.  The 
range  vector  is 

R_'  -  i/  ( |  R_|  cos  p  sin  p) 

+  j/  f(|R|  cos  p  cos  p)cos  6  —  (|  R|  sinp  sin  6] 

+  k7  f(|  R  |  cos  pcos  p)sin  6  +  ( |  R_ |  sin  p)cos  5]  (5) 


and 
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which  expression  reduces  to  Eq.  4  for  5  =  0.  The  new  expressions  for  the 
pitch  and  roll  angles  are 


sm  p  = 


R./ 

1 


and 


J 2  2 

V  R./  +  R ./ 

1  J 


R 


tan  p  = 


ky 


Vr.,2  +  R./2 

1  J 


The  radius  vector  r_  to  the  piercing  point,  P  is 


where  I  r  I  =  6378.  16  kilometers.  The  coordinates  of  p  are 


tan  latitude  = 


K 


2  2? 


and 


sin  0 


Vrf  +  rj 


J- 2~  ^ 

Vri  +rj 


(6) 


(7) 


£  =  J^(|RL|cosp/sinp/) 

—  (m|j;|  cos  6  —  |R|  cos  p  cos  p7) 

+  K  ( I  R  I  sin  p7  —  m  I  r  I  sin  6)  (8) 


(9) 


(10) 


The  longitude  (East)  of  P  follows  from  Eq.  10  and  from  the  longitude  (East), 
0,  of  the  sub-satellite  point: 


Longitude  (East)  =  0  +  07 


(11) 


In  Fig.  1  the  coordinates  have  been  arbitrarily  aligned  such  that  the  sub-satellite 
longitude  is  270°E  . 
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Having  computed  the  latitude  and  longitude  of  the  piercing  point  of  II, 
one  would  like  to  form  an  arbitrary  number  of  vectors  all  lying  on  the  surface 
of  the  cone  and  similarly  compute  their  piercing  points.  Let  the  coordinate 
system  Ox/y/z/  be  rotated  by  the  angle  p/  around  the  axis  OZ '  ,  then  the 
resultant  system  be  further  rotated  by  p  around  the  resultant  x  axis.  The 
transformation  is: 
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In  the  doubly-primed  system  let  there  be  a  vector  u  as  illustrated  by  Fig.  2. 
The  magnitude  of  u  is 

|  u  |  =  |  Rj  tan  a  (13) 

and  the  instantaneous  value  of  u  is 

u  =  i"  |u|  cos  (3  +  k"  |u|  sin  (3  (14) 

where  (3  is  a  convenient  angular  parameter,  measured  clockwise  from  the 

f  Vi 

x"y"  -plane.  The  n—  M  vector,  M  ,  lying  on  the  surface  of  the  cone  is 
defined  as 


M  =  R  +  u 
— n  —  — n 


(15) 


which  can  be  written  in  the  single -primed  coordinate  system  as 
- 2 - Z - 21 

R./  +  R./  +  R,  /|i/(cos  p/sinp/  +  tan  a  cos  (3  cos  p7-tan  a  sin  P  sin  p* 
ijk[—  v  r  n  n 

+  ^(cos  p'  cos  p f  —  tan  a  cos  sin  p7  —  tan  a  sin  (3^  sin  p/  cos  p7) 

4  k/(sin  p  V  tan  a  sin  (3  cos  pOj  ;  (16) 


sin  p') 
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and  in  the  Oxyz  coordinate  system  as 


/  2 2  2*  ( 

M  =wR.  /  +  R.  /  +  R,  /  |_i(cos  p 7  sin  p7  +  tan  a  cos  P  cos  p7  —  tan  a  sin  (3  sin  p7  sin  p7) 
n  *  ^  j  k  1  n  n 

+  j_  [(cos  p/  cos  p7  —  tan  a  cos  (3^  sin  p7  —  tan  a  sin  (3^  sin  p7  cos  p7)  cos  6 

+  (sin  p7  +  tan  a  sin  (3  cos  p7)sin  6 

n  J 

+  kj(cosp/cosp/  —  tan  a  cos  P^  sin  p7  —  tan  a  sin  P^  sin  p7  cos  p/)(  —  sin  6) 

s]  J  .  (1 7) 


+  (sin  p7  +  tan  a  sin  P  cos  p7)  cos  S| 


From  Fig.  2  the  following  trigonometric  relations  can  be  written: 
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cos  b 
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1  2  2~ 1 

M  .  +  M  .  +  M  . 
m  nj  nk 


(18) 
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m  sin  b 

(19) 
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have  two  values:  c ^ 

o 

o 

O' 
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cn2  =  180°  -  cnl 

>  90°  ; 

normally,  c ^  should  be  used  only  (unless  one  is  interested  in  the  inter¬ 
section  of  the  sphere  and  the  cone  emerging  from  the  sphere).  The  expression 
for  pitch  and  roll  angles  associated  with  M ^  are 
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The  radius  vector 


r 

— n 


to  the  piercing  point  of  M  ^  is 


r  =  I(M  cos  p  sin  p  ) 

— n  — v  n  n  *n 

—  J(m|r|  cos  6  —  M  cos  p  cos  p  ) 

—  i—1  n  n  n 

+  K(Mn  sin  Pn  “  mll.l  sin  6)  (24) 


where  M  is  obtained  from  Eq.  21,  M  /  I M  I  of  Eq.  17.  The  coordin- 
n  n  1  — n 1 

ates  of  P  are 
n 


tan  Lat 


nK 


n  V 


(25) 


2  2 

r  T  +  r  T 

nl  nj 


sin  0 


'nl 


n 


V 


- 2 - T 

r  r  +  r 
nl  nj 


(26) 


Longitude  (East)  =  0  +  0' 
n 


(27) 


Equations  (25)  and  (26)  yield  the  latitude  and  longitude  coordinates  of 

the  point  where  the  vector  M  ,  which  is  at  an  angle  a  to  the  vector  R 

n  — n  °  — 

and  has  the  same  origin,  pierces  the  surface  of  the  sphere.  The  M  vectors 

—  n 

all  lie  on  the  surface  of  the  cone,  their  piercing  points  on  the  sphere  are  all 
points  on  the  curve  of  first  intersection.  By  suitably  choosing  values  of  (3^  , 
the  latitude  and  longitude  coordinates  of  any  point  on  this  curve  C  can  be 
readily  computed;  hence,  C  can  be  described  with  any  desired  degree  of 
precision.  Many  antenna  power  patterns  are  circularly  symmetrical,  this 
solution  to  the  cone-sphere  intersection  problem  is  directly  applicable  in 
these  instances.  The  method  presented  here  could  be  readily  extended  to 
other  antenna  pattern  shapes  and  cross-sections  such  as  ellipses.  All  that 
is  needed  is  the  replacement  of  u  --which  is  but  a  parametric  equation  of 
the  circle  (a  conic  section)  with  (3  as  the  parameter  - -with  a  suitable  para¬ 
metric  expression  describing  the  shape  of  the  cross-section  of  the  antenna 
pattern. 
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A  practical  variation  of  the  problem  presented  above  is  the  situation 
when  the  range  vector  is  unknown  but  the  latitude  and  longitude  of  the  spot  P 
is  given  toward  which  the  antenna  boresight  is  directed.  In  that  case  Eq.  4 
is  replaced  with 

R_  =  _i(|r_|  cos  Lat  sin  07) 

+  j_  |(m|rj  cos  ^  —  |  rj  cos  Lat  cos  07)  cos  6 
+  (m|r_|  sin  6  +  |r_|  sin  Lat)  sin  fij 

+  k  j(m  |  rj  cos  6  —  |rj  cos  Lat  cos  07)  (  -  sin  6 ) 

+  (m|rj  sin  6  +  |rj  sin  Lat)  cos  ,  (28) 

and  Eq.  5  is  replaced  with 

R*  =  ±'(\l\  cos  si-n  07) 

+  j^mlr^l  cos  6  —  |r_|  cos  Lat  cos  <fi) 

+  ^(mlr^l  sin  6  +  |r_|  sin  Lat)  .  (29) 

The  solution  proceeds  as  above  with  the  aid  of  Eqs.  6  through  27.  The  first 
results,  the  coordinates  of  the  piercing  point  calculated  with  the  aid  of 
Eqs.  9-11  will  produce  the  same  values  which  were  given  for  the  coordinates 
of  the  spot  P  .  The  subsequent  equations  will  yield  coordinates  of  points  P^ 
on  the  curve  of  intersection  C  . 

In  the  preceeding  discussion  it  was  tacitly  assumed  that  the  first  inter¬ 
section  was  on  the  surface  of  the  sphere  in  its  entirety.  Let  us  next  consider 
the  question  of  coverage.  The  subtended  angle,  ^max>  the  Earth  from 
the  origin  of  Oxyz  is 

sin  (b  )  =  -  (30) 

max  m 


9 


(31) 


One  of  the  direction  cosines  of  the  range  vector  is 


cos  b 


R. 

1 


There  are  three  possibilities:  (a)  Piercing  point  of  boresight  is  on  surface 

of  the  Earth,  at  least  one-half  of  antenna  beam  illuminates  the  Earth 

fb  <  b  ]  ;  (b)  Range  vector  is  tangent  to  surface  of  the  Earth,  at  most 

max 

one -half  of  antenna  beam  illuminates  the  Earth  [b  =  ^max]  '>  (c)  Piercing 

point  of  boresight  is  not  on  surface  of  Earth,  limited  coverage  is  still  pos¬ 
sible  [b  >  t>max]  •  These  statements  are  valid  under  the  assumption  that 
the  vertex  angle  of  the  cone  is  about  the  same  as  the  angle  subtended  by  the 

Earth  as  seen  from  the  position  of  the  satellite.  In  case  of  the  M  vectors 

— n 

on  the  surface  of  the  cone,  the  case  of  interest  here  is  when  b  >  b  for 

n  max 

some  value  of  b^  .  In  that  case,  one  is  interested  in  the  latitude  and  longitude 
coordinates  of  the  points  on  the  curve  of  last  contact  between  the  sphere  and 
some  inside  portion  of  the  cone.  The  computations  proceed  as  follows. 

Replace  Eq.  24  with 


r 

— n 


r  J  I(cos  b  sin  t  ) 
— 1  )—  max  n 


+  J(  —  cos  6  sin  b 

—  max 


+  K(  —  sin  6  sin  b  + 

—  max 


sin  6  cos  b  cos  t  ) 
max  n 

cos  6  cos  b  cos  t  )i  (32) 
max  n  / 


and  continue  with  Eqs.  25  through  27  to  find  the  coordinates  of  the  points  on 

the  curve  of  last  contact.  In  Eq.  32  the  tilt  angle,  t,  defined  positive  when 

measured  counter-clockwise  from  yz-plane  (towards  yx-plane),  between  the 

yz-plane  and  the  plane  of  the  triangle  formed  by  the  vectors  r  ,  mr  ,  M 

— n  —  — n 

is 


cos  t 
n 


A±=n\ 

yli  x  M.n  y 


(33) 
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This  completes  considerations  of  the  general  problem  of  intersection  of  a  cone 

and  a  sphere  with  the  exception  of  the  trivial  case  when  b  >  b  for  all 

n  max 

values  of  b  .  In  this  last  case  there  is  no  intersection  and  either  the  entire 
n 

hemisphere  is  illuminated  or  no  portion  of  the  surface  of  the  Earth  is  illum¬ 
inated. 


III.  EXAMPLES 

As  an  illustration  of  the  preceeding,  consider  the  following  examples. 


Subsatellite  Longitude,  0 
Cone  vertex  angles,  a 
Increments  of  angular  parameter  (3 
Coordinates  of  piercing  point  of  boresight 


270°E 


0.6,  1.0,  1.2 

,o 


10' 


288. 733  E 
42. 462°N 


Distance  to  satellite  from  center  of  Earth, m^  6.611  Earth  radii 


Instantaneous  declination,  6 


2.0 


The  results  of  the  computation  are  shown  graphically  in  Fig.  3.  In  the 
second  example,  shown  in  Fig.  4,  the  same  parameters  are  used  as  above 
except  mr  =  19.  832  Earth  radii,  and  (3  =  1°  . 
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Fig.  3.  Footprint  of  a  narrow  beam  antenna  from  synchronous  orbit. 
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Fig.  4.  Footprint  of  a  narrow  beam  antenna 
from  a  distance  of  19  Earth  radii. 
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